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I— I Presently available perturbative QCD calculations combining hard process matrix element with 

the Parton Shower Monte Carlo programs feature hard process matrix element calculated often 
beyond the leading order (LO), that is including complete next-to-leading-order (NLO), or even 
next-to-next-to-leading-order (NNLO) corrections, while Parton Shower is only at the leading 

C ■ order (LO). We report here on a work in progress which demonstrate feasibility of constructing 

Parton Shower Monte Carlo (PSMC) featuring complete NLO corrections to QCD evolution with 
respect to the logarithm of the factorization scale. This effort presently covers non-singlet subset 
of Feynman diagrams contributing to the above QCD evolution. It should be stressed that our 
. • approach to the NLO QCD evolution is exclusive, that is giving insight into the fully unintegrated 

. ^ phase space. However, at the inclusive level our implementation agrees exactly with the standard 

inclusive picture of the NLO DGLAP evolution. Our new approach (after including NLO singlet 
diagrams) provides a complete method of combining the resummed and fixed order perturbative 
QCD calculations beyond the LO in a form suitable for the MC implementation, an alternative 
to the existing ones. First practical applications will include Monte Carlo generators for W/Z 
production processes in hadron-hadron colliders and lepton-hadron colliders. 
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As the Large Hadron Collider (LHC) just started to operate, experiments will need high ac- 
curacy QCD calculations not only for inclusive observables, but also for multiparton final states, 
especially in a form of events from Parton Shower Monte Carlo (PSMC) programs. Huge range of 
energy scales in the PSMCs will require including beyond the LO calculations not only in the hard 
process part but also in the middle of the shower (ladder). Why 25 years old traditional PSMCs, see 
refs. [ 1 ] , are still restricted to the LO level only? Main reasons are that: (A) principal role of PSMCs 
was to hadronize quark and gluons into hadrons, (B) hadron collider data were until recently rather 
poor, (C) conceptual problem on the side of the factorization theorems of QCD [2], badly suited 
for the MC implementation, as they (i) violate 4-momentum conservation, (ii) irreversibly inte- 
grate over real emissions (iii) feature huge over-subtractions and cancellations (iv) operate with 
non-positive distributions, (D) last not least computers were until recently to slow to consider unin- 
tegrated NLO corrections inside the PSMC ladder. With the presented works an era of NLO level 
PSMCs finally comes soon! Possible profits/gains from NLO PSMC are the following: (a) new 
simpler and cleaner schemes of matching of the hard process ME at NLO and NNLO level with 
NLO PSMC, (b) natural extensions towards BFKL/CCFM evolution at low x, (c) better modelling 
of low scale phenomena, Q < lOGeV, (d) new MC tools for porting information on parton distri- 
butions from one process to another (from DIS at HERA to W/Z production at LHC) and from one 
scale to another - an alternative for traditional inclusive PDFs. In our work presented here, we take 
scheme of Curci, Furmanski, Petronzio (CFP) [4] as a guide and we concentrate on the non-singlet 
NLO evolution for the initial state radiation (ISR) of multiple gluons. NLO corrections are added 
to each LO kernel in the middle of LO ladder of the PSMC, see below. 



XT' 




such that at the inclusive level we reproduce NLO DGLAP [3] evolution. As depicted above, our 
starting point is the LO ISR single ladder like in DIS. Let us start with including NLO conections 
to the gluon emission vertex at the top of this LO ladder. 
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where the Monte Carlo weights reweighting LO distributions to NLO level is given by: 



1. 



(2) 



Let us explain first all parts and the ingredients in Eq. (1). First term (ladder) is just the LO ladder 
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multiplied by NLO virtual correction factor for the last vertex, where {kij — — — 
LO Altarelli-Parisi evolution kernel, QQ>ai>ai-i indicates the ordering in angle (a, is the rapidity 
related variable defined as a,- = |k,j^|/a,-, a,- is plus variable in Sudakov parametrization), ^ is the 
Lorentz invariant phase-space and e^^'^x is the Sudakov formf actor. The NLO virtual coiTcction 
■■■■ = (l +29t(Aj^j)) t simply multiply the LO Born (last) vertex and since it is finite, it can 
DC easily implemented by reweighting the LO distribution. The real corrections coming from the 
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with the soft counterterm 



subtracted bremsstrahlung type diagrams j |[|' | = j J , + j ~ j ' 

I subtracting the LO component affect pairs of partons (act in two particle phase-space) ^ The 
distinctive feature of Eq. (1) is the presence of the summation over trailing gluons along the ladder. 
Without this summation (taking only W {kn,kn-\)) we wouldn't have the full NLO correction. 

Including the full NLO correction all over the ladder requires to sum over the positions of the 
real NLO insertions (corrections) in the ladder and adding virtual corrections to all vertices: 
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The above scheme is tested numerically with 3-digit precision using prototype MC program: 




I NLO INCLUSIVE | 




■0.5 0.5 1 1.5 2 2.5 3 3.5 4 4.5 -0.5 0.5 1 1.5 2 2.5 3 3.5 4 4.5 



' In our graphical notation the amplitude square means cut diagrams for both Feynman graph and the counterterm. 
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Perfect agreement of the new scheme and traditional inclusive approach is seen, especially in the 
left lower plot, where the ratios for two sums in Eq. (3) are plotted separately. In the right plot we 
see the weight distributions for inclusive and exclusive approach. The exclusive weight is positive 
and naiTow. 

In NLO corrections proportional to CaCf color factor, we get pair production diagram | . 
First question is why we treat this type of corrections in a special way not just include them in the 
same manner as before. The problem is that it induces Sudakov double logarithm +Sfsr (addi- 
tional singularity when produced pair is coUinear) which would ruin the MC weight, if we would 
restrict ourself to strict NLO only. Resummation/exponentiation for this contribution is mandatory, 
and new LO ladder on the external leg has to be implemented in the MC - each emitted gluon will 
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cascade into gluon emission ladder. The FSR counterterm f is constructed, encapsulating 
exactly soft and coUinear singularities of C^Cf diagrams (both double and single logs). The ex- 
pression for the MC distribution with the last emission in the ladder upgraded from the LO to NLO 
level reads: 
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where the MC weights are given by 
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are finite, since all coUinear and soft, ISR and FSR, singularities are subtracted. In particular, in the 

II I' i+m' 

virtual correction = ( 1 + 29t (A^^,^ + V^^^ — 5^^^)) i-'-- FSR Sudakov 5^^,^ is subtracted. 

Three-digit precision numerical test of the above FSR methodology was done for a single NLO 
ISR-i-FSR insertion for n = 1 , 2 ISR gluons and infinite number of FSR gluons. 
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The comparison is done between the MC and analytical results. 

The exclusive approach presented here has been tested against an independent MC implement- 
ing standard inclusive DGLAP evolution in the NLO level. We can see the perfect agreement of 
both approaches on the level of few promiles, see [5, 6, 7]. Let us stress that in our exclusive MC 
scenario we have complete insight into the transverse degrees of freedom up to NLO level, some- 
thing that is unavailable in other common inclusive or semi-inclusive approaches. Our exclusive 
NLO MC features good behaviour of the MC weights - weights are positive and narrowly dis- 
tributed. The techniques presented above are promising and will be exploited in the construction of 
NLO MC generators for W/Z production @LHC and DIS @HERA. Works related to this ongoing 
project are reported in refs. [5, 7, 8, 9, 10]. 
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